Dissipationless mechanism of skyrmion Hall current in double-exchange ferromagnets 
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We revisit a theory of skyrmion transport in ferromagnets. On a basis of an effective U(l) gauge 
theory for spin-chirality fluctuations in double-exchange ferromagnets, we derive an expression for 
the velocity of a skyrmion core driven by the dc electric field. We find that mutual feedback effects 
between conduction electrons and localized spins give rise to Chern-Simons terms, suggesting a 
dissipationless mechanism for the skyrmion Hall current. A conventional description of the current- 
induced skyrmion motion, appearing through the spin transfer torque and scattering events, is 
reproduced in a certain limit of our description, where the Chern-Simons terms are not fully incor- 
porated. Our theory is applicable to not only metallic but also insulating systems, where the purely 
topological and dissipationless skyrmion Hall current can be induced in the presence of an energy 
gap. 

PACS numbers: 



I. INTRODUCTION 

Soliton dynamics in the presence of fermions has been 
one of the most fundamental issues in various fields of 
physics. It plays a central role in domain-wall dynamics 
in conducting polymers [l| and vortex dynamics in super- 
conductors in the field of condensed-matter physics [2| . It 
is also relevant to a confinement and baryon dynamics in 
high-energy physics Q . Skyrmions [3| were shown to ap- 
pear as elementary excitations in the quantum Hall sys- 
tem and the interplay between skyrmions and spin- 
wave excitations was also studied @ . Actually, skyrmion 
excitations have been observed in the cold atom system 
0, and their crystallization has also been observed in 
both two and three dimensional helical magnets [1, Q . 

Advances in spintronics have promoted intensive and 
extensive studies on domain-wall dynamics in magnetic 
systems for its application to the magnetic memory de- 
vice [T]| . In particular, spin-polarized electric cur- 
rents allow for an efficient control of domain- wall dynam- 
ics because of the spin torque transfer. A conventional 
theoretical approach to this phenomenon is based on the 
Landau-Lifshitz-Gilbert equation, that is, the equation of 
motion for a single spin. In this formalism the so-called 
Gilbert damping term is introduced in a phenomenolog- 
ical manner jllf to take account of the dissipation in the 
spin dynamics, caused by the coupling to the conduction 
electrons and the relativistic spin-orbit coupling. It re- 
veals that the spin current in the itinerant ferromagnet 
drives a domain- wall motion to the longitudinal direction 
[ll| . The spin dynamics in nanoscale magnets has also 
been studied by introducing collective coordinates [l^ . 
such as positions of domain wall and magentic vortices. 
Then, it has been argued that the spin current generates 
half magnetic vortices and/or anti- vortices which can ex- 
hibit a nontrivial motion, including the transverse motion 
in the presence of the Gilbert damping. In these previous 



approaches spin currents are given, thus feedback effects 
of the spin dynamics on the electron dynamics have not 
been considered. 

In this paper, we develop a theory of skyrmion trans- 
port in double-exchange ferromagnets, where conduction 
electrons interact with localized spins via the Hund's- 
rule coupling. An essential aspect of our study lies in 
a mutual feedback effect between electron spin currents 
and skyrmion currents, resulting in the dissipationless 
skyrmion Hall current of the topological origin, which 
should be distinguished from the dissipative skyrmion 
Hall current replying on the scattering events [ll| . Actu- 
ally, a coupling between the spin current and a skyrmion 
motion has also been discussed in Refs. [l3l - [l6j . where 
topological magnetic textures generate electric currents 
via the Berry-phase induced electro-motive force, which 
produces feedback on the magnetization dynamics via 
a spin-transfer torque. In these pioneering works, how- 
ever, they considered the dissipative mechanism of the 
skyrmion Hall current but not the dissipationless mech- 
anism. 

In Sec. ini resorting to the CP^ representation for 
the localized spin, we derive an effective field theory 
for itinerant electrons and bosonic spinous, which inter- 
act via gauge fluctuations representing the spin chirality. 
Based on this effective field theory for the strong Hund's- 
rule coupling limit [ll| . we derive the Maxwell-Chern- 
Simons equations for both internal U(l) gauge fiuctua- 
tions and external electromagnetic fields. The emergence 
of Chern-Simons terms is ascribed to mutual feedbacks 
between itinerant electrons and skyrmions. Focusing on 
the center-of-mass motion of the skyrmion, we obtain the 
velocity of the skyrmion core in terms of the dc electric 
field. The Chern-Simons terms induce the dissipation- 
less skyrmion Hall current normal to the applied electric 
field. For comparison with previous theories, we show 
m Sec. HH] that our U(l) gauge-theory formulation for 
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skyrmion dynamics reproduces the dissipative skyrmion 
Hall current in the Landau-Lifshitz-Gilbert-equation ap- 
proach. 

In fact, the topological contribution of the dissipation- 
less skyrmion Hall current scales with the skyrmion den- 
sity, and thus vanishes in the thermodynamic limit when 
a single skyrmion is considered. However, we find that 
the Rashba spin-orbit coupling for conduction electrons 
produces a finite skyrmion Hall current in the thermody- 
namic limit, allowing for an observation of this intrinsic 
skyrmion Hall effect. This could be realized at the sur- 
face of three-dimensional topological insulators when 
a ferromagnet is deposited. Recently, a dissipationless 
mechanism for magnetization switching was proposed in 
the topological surface state, where the Chern-Simons 



term plays an essential role \L 



II. U(l) GAUGE THEORY FOR SINGLE 
SKYRMION DYNAMICS 

A. An effective Maxwell-Chern-Simons Lagrangian 

We start from an effective U(l) gauge-field formulation 
of a two-dimensional double exchange model with the 
Rashba spin-orbit coupling, where itinerant electrons in- 
teract with localized spins via the ferromagnetic Hund's- 
rule coupling Jh- This is described with the following 
partition function Z and the Lagrangian density £; 



Cb = 2iSar, = i'lidr - ih- - JhScf - iaar - iAr)'ijja + t\[dr - ic^a^ - iAr - ia{Xso/t)zla'^pZi3]ipo^\'^ , 

Cz = Pszlidr - iar)Za ^ tps\[dr - Zttr ^ i{\so / t)z'lal^pZ l3\Za\'^ . (1) 
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This effective field theory can be derived along the con- 
cept of Ref. The derivation is given in Ref. [l^l 
as well as in Appendix A. The underlying mechanism to 
justify this formulation is that the spin dynamics of itin- 
erant electrons instantaneously follows that of localized 
spins in the strong Hund's-rule coupling limit as far as 
the dynamics of localized spins are much slower than that 
of itinerant electrons fllj . Below we explain the physical 
meaning of each term. 

Lb represents the single-spin Berry phase, resulting 
from the curved nature of the SU(2) spin manifold. This 
term is indispensable for reproducing the skyrmion dy- 
namics that has been obtained in the Landau-Lifshitz- 
Gilbert equation approach flll |. 

describes the dynamics of itinerant electrons, where 
ipa and Za represent the fermionic field for itinerant elec- 
trons and the CP^ spinon field for localized spins, respec- 
tively. As mentioned above, spins of itinerant electrons 
follow those of localized electrons in the strong Hund's- 
rule coupling limit. This constraint provides dynamics 
of itinerant electrons with an effective internal fiux, orig- 
inating from the curvature in the spin space. Namely, 
their orbital motion is affected by an effective Aharonov- 
Bohm phase or the Berry-phase connection, described by 
an internal U(l) gauge field with /i = T,x,y. Phys- 
ically, this gauge field represents spin chirality fiuctua- 
tions, and couples to electrons with opposite signs of cou- 
pling constants for "spin" up and down. In contrast, 
represents the external electromagnetic potential which 
linearly couples to electric charge/current density of itin- 
erant electrons, p.^ is the chemical potential, which is 
determined to fix the total number of itinerant electrons. 



t is the hopping energy of itinerant electrons. The last 
term represents a spin vector potential, which originates 
from the Rashba spin-orbit coupling Ago- This provides 
an interaction between the spin and the "spin current" , 
and thus quenches the spin direction to that of the mo- 
mentum or spin current. 

Cz describes the dynamics of the CP^ spinon field Za 
for localized spins, in other words, their directional (an- 
gular) fiuctuations. In ferromagnets, this spinon dynam- 
ics produces the magnon excitations, which exhibits the 
u! (X k'^ dispersion relation in the SU(2) symmetric case. 
We have introduced the spin density of itinerant elec- 
trons, ps = l^Ylia'^'^t'i''^'^ ■ We have an additional term 
generated by the Rashba spin-orbit coupling Ago- 

Note that the skyrmion configuration creates a non- 
trivial background potential for the Berry gauge connec- 
tion and thus a fictitious internal magnetic field in 
the z direction. This affects the dynamics/transport of 
itinerant electrons. Indeed, even if the skyrmion is static, 
it produces the anomalous Hall current of itinerant elec- 
trons [21I, . Actually, this is one side of the mutual 
feedback effect between the skyrmion and the fermionic 
matter. On the other side, the topologically induced 
anomalous Hall current is accompanied by the dissipa- 
tionless skyrmion Hall current, when the skyrmions are 
depinned intrinsic objects, as we will show later. 

The skyrmion motion can be uncovered from an ef- 
fective action for both U(l) Berry gauge fields and 
electromagnetic fields A^, integrating over electrons 
and spinons z„ in the skyrmion background. We separate 
the Berry gauge field into two pieces which correspond to 
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its classical configuration and quantum-fluctuation part, 
respectively. The classical configuration of the Berry 
gauge field is determined from an equation of motion 
for spinous, where dynamics of spinous is taken into ac- 
count classically. The skyrmion solution of the spinon 
field gives rise to an effective magnetic field for electrons, 
given by the following relation of = —^[zp{dfj,z^) — 
{df^zp)z^], where the superscript c denotes "classical". 
Inserting the effective magnetic field into the Schrodinger 



where i5a^ represents the quantum-fluctuation part. 

In the Maxwell Lagrangian Cm, o'ss, cTcc, <^sc, and 
Xss, Xcc, Xsc are conductivities and diamagnetic sus- 
ceptibilities, associated with spin-current-spin-current, 
chargc-currcnt-charge-current, and spin-current-charge- 
current correlation functions, respectively. Electrons are 
assumed to be in the diffusive regime, resulting in the 
z = 2 dynamics for gauge fluctuations, where z is the 
dynamical exponent to represent the dispersion relation, 
uj (X . is the projection operator for transverseness 
of the gauge dynamics, given by P^j = Sij + didj/{—d^), 
where i,j = x,y and —d"^ = —d^ — dy. Dynamics of 
the temporal part Sa-r can be neglected in the low energy 
limit because such fluctuations are gapped and decoupled 
with spatial fluctuations in the Coulomb gauge. 

In the Chern-Simons action Ccs, 0ss(x — X) denotes 
a local spinon-Hall conductance, given by the transverse 
spin-current-spin-currcnt correlation function. Namely, 
it describes the spin Hall current generated by the mag- 
netic field gradient. 0cc(x— X) is a local charge- Hall con- 
ductance, given by the transverse charge-current-charge- 
current correlation function. 0sc(x — X) is a local spin- 
Hall conductance, given by the transverse spin-current- 
charge-current correlation function. X represents the 
spatial coordinate of the skyrmion core, which should be 
distinguished from that of the fields, i.e., x. Actually, the 
skyrmion contribution to the charge Hall current decays 
with the distance |x — X| from the skyrmion core. In 
particular, the charge Hall current vanishes at the long 
distance if the relativistic spin-orbit coupling Ago is ab- 
sent. All coefficients in this effective field theory will be 
found in section IV. 

Based on Eq. ([2]), we investigate the skyrmion dynam- 
ics under an external dc electric field applied along the x 
direction. Since we are interested only in the transport 
properties in the linear response to the applied dc clcc- 



equation for electrons, we construct the space of wave 
functions, well fitted to the skyrmion potential. Then, 
we can integrate over electrons, and expand the result- 
ing logarithmic term up to the second order for gauge 
fluctuations. The whole procedure is shown in section 
IV. This gives rise to not only the Maxwell Lagrangian 
Cm but also the spatially dependent Chcrn-Simons ac- 
tion Ccs for the gauge-field dynamics. 



(2) 



trie field, it is sufficient to treat the constant velocity of 
the skyrmion core. To facilitate the calculation on the 
coupling between itinerant electrons and localized spins, 
it is convenient to introduce the frame moving with the 
skyrmion core at the origin. In this skyrmion moving 
frame, the time derivative and the time component of 
the gauge field are transformed as 

dr > dr — ■ Or, > Sa^- — • Ssir, (3) 

where is the constant skyrmion velocity driven by 
the external electric field, which will be determined self- 
consistently below. Note that we implicitly ignore the 
modification of the shape of the static single-skyrmion 
configuration, which requires more careful self-consistent 
treatment of the coupling between itinerant electrons and 
localized spins, but does not spoil the topological origin 
of our skyrmion Hall current totally. 



B. Skyrmion dynamics under electric field 

Taking the derivative of Cef / with respect to A^^ and 
6a^, we obtain an equation of motion for the U(l) Berry 
gauge field and that for the electromagnetic field, respec- 
tively, 

Xsd^Sa^ + Xcsd^A, = 2(3 - M)v, 

+ assSci —eijSej - -^e^jEj, 

+ accEi —f-ijEj - — ^e.y(5ej(4) 

TT ZTT 

where the imaginary time has been replaced with the real 
time. The terms linearly proportional to the skyrmion 
velocity Vi originate from the Berry-phase term in the 
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moving frame. M ^ criplipa + PsJ^a ^l^<^ 

responds to the magnetization density, which effectively 
reduces the coefficient of the Berry phase term and sup- 
presses that of the skyrmion velocity, pe/ ^JV'cr^ 
is the charge density. Ej is an external electric field and 
6ej = ejfj,,/dfj,Sai, is an internal electric field. 

Physics of these Maxwell-Chern-Simons equations can 
be understood as follows. Recalling the structure of 
the Maxwell equation, one can construct two constituent 
equations, which relate "spin" and "charge" currents 
with both internal and external electric fields, 



2tt 

^Jfa = ^ccE^ - ^^^]EJ - ^eijSCj, (5) 



where jf^ represents the current of itinerant electrons 
with the spin index a flowing in the j-direction, given by 



(6) 



The spin conductivity ass vanishes in the paramagnetic 
phase without the spin-orbit coupling, while in the fer- 
romagnetic phase with the Zeeman splitting JhS, it is 
finite and the spin current of itinerant electrons is gen- 
erated by the external electric field Ei. Equation ^ 
generalizes the standard constituent relation in metals, 
where dynamics of conduction electrons are in the diffu- 
sive regime, introducing the Chern-Simons contribution 
into the equation, which plays an essential role for the 
mutual feedback effect between skyrmions and itinerant 
electrons. The presence of the Chern-Simons term in the 
constituent equation confirms the anomalous Hall effect 



of itinerant electrons in the metallic ferromagnet, as dis- 
cussed before. 

It is straightforward to solve coupled equations 
when the spatial dependence for both gauge fields is ne- 
glected. Performing integration of J dx J dy in Eq. Q , 
we obtain 



-2{S - M)v., = assSe., - ^U,5e, ~ ^e^jE,, 
a" a" 
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(7) 



where 



J '^^ J '^yQ'^-'sW' 
'^cc= J dx J dyQcci^), 

CT^ = J dx j dyQsci'^) 



(8) 



are Hall conductivities with L being the linear spatial 
dimension of the system. 

Coupled equations (O describe both the internal elec- 
tric field 5ei and the skyrmion velocity Vi as a function 
of the external electric field Ei. We find the following 
expression for the internal electric field 



Pel 
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Inserting this expression into Eq. we obtain the 

skyrmion velocity as a function of the external electric 
field E^, 



{S - M)ag + peia^s) ((^ss<y^c + ^ccT,2 I + PeKJss U^L^CJssCJcc - (J^s^. 



(< 



Ex 



((5 - M)af, + peia^^y + n^L^p^^ass)' 



E..r. 



(10) 



To understand the above expression, we consider two 
limiting cases. First, we take the limit of crf^ = cr^ = 



0; resulting in 

Or 



-Ex 



-Ex 
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(11) 



The skyrmion current is driven to not only the same di- 
rection as the applied electric field but also the orthogo- 



nal direction corresponding to the Hall motion. 

Second, we take another limit of Uss = fee = cr.sc = Oi 
corresponding to an insulator. Then, we find 



0, 



Ex 



(12) 
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This is a remarkable result. Although the skyrmion Hall 
current induced by the electric field vanishes in the ther- 
modynamic limit, the nature of the skyrmion Hall cur- 
rent is dissipationless. We cannot find any coefficients 
associated with conductivity, giving rise to dissipation. 
This certainly originates from the Chern-Simons terms, 
an essential feature of the "self-consistent" treatment. 

It is interesting to observe that the dissipationless 
skyrmion Hall current may not vanish in the thermo- 
dynamic limit if the spin-orbit interaction is introduced. 
The Hall coefficient without the spin-orbit coupling is 
proportional to the density of skyrmions, thus ~ in 
the case of the single skyrmion. As a result. Hall conduc- 
tivities of afg, CT^, and crf^ in Eq. (fT2l) are constants even 
after the spatial integration. However, the spin-orbit in- 
teraction can give rise to a finite value for the Chern- 
Simons coefficients in the thermodynamic limit, thus its 
integral value corresponding to cr^, cr^, and af^ will be 
proportional to L^. Then, we obtain the dissipationless 
skyrmion Hall current in ferromagnetic insulators with 
the spin-orbit coupling. 



C. Discussion : Dissipationless skyrmion Hall 
current in the surface state of three dimensional 
topological insulators 

We discuss the dissipationless skyrmion Hall current 
in the surface state of three dimensional topological in- 
sulators, where magnetic impurities are deposited. An 
effective field theory for surface Dirac electrons with lo- 
calized spins is given by 



current 



1 



Jn = —ef^Luxfi ■ dyfi X d\n 



Sir 



corresponds to the skyrmion current, discussed before. 
The minimal coupling between the skyrmion current and 
the electromagnetic field implies that a skyrmion car- 
ries an electric charge, where a normalizable fermion 
zero mode exists on the topological soliton, inducing a 
fermionic charge. On the other hand, the last term shows 
the geometric phase, identified with the Hopf term, which 
determines the statistics of the topological soliton and its 
spin quantum number. These two terms are well known 
in the field theory, referred as the quantum anomaly, 
where the induced fermionic charge for a soliton is a phe- 
nomenon due to the local anomaly while the topological 
phase is due to the global anomaly [1^ [IJ] ■ Mathemat- 
ically speaking, the existence of such topological terms 
is guaranteed by the fundamental property of the elliptic 
operator, the Dirac operator in the present case, called 



the Atiya-Singer index theorem [25 



One can express the above effective field theory as fol- 
lows, resorting to the CP^ representation. 



Seff = Se 



X 



An 



2e2 



(15) 



with 

Ju 



S = Se 



'1- 



Sb is the single-spin Berry phase term to appear from 
the coherent-state representation for the localized spin n 
in the path-integral quantization. -0 represents the sur- 
face Dirac fermion in the irreducible representation, and 
the covariant derivative is iD = 7^(i9^ -I- A^)^ where 7^ 
with /i = r, a;, y is two by two Dirac matrices and is an 
external electromagnetic vector potential, m > is an ef- 
fective coupling constant between surface Dirac fermions 
and deposited magnetic impurities. Compared with the 
double exchange model, the only difference is that non- 
rclativistic electrons are replaced with Dirac fermions. 

Integrating over Dirac electrons and performing the 
gradient expansion for the resulting logarithmic term 
[23| . one finds the following effective action 



Seff = Se 



1 



iA^J^ 



(14) 
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where dynamics of localized spins is governed by the non- 
linear (T model with additional terms. The conserved 



g cx 1/m, where the skyrmion current is 
f_i — -^^fiuxdvO-x and the topological phase is 'd[fi] = 
j^^fii^xdfj.duO'X- This effective field theory for dynamics 
of localized spins on the topological surface is essentially 
the same as our previous effective field theory for dynam- 
acs of localized spins in the presence of non-relativistic 
electrons except for the fact that coefficients in topolog- 
ical terms, i.e., Chern-Simons terms have finite values in 
the thermodynamic limit. As a result, we find the dissi- 
pationless skyrmion Hall current in the surface state of 
three dimensional topological insulators, although sur- 
face Dirac electrons become gapped and insulating due 
to time reversal symmetry breaking. Furthermore, the 
electric charge of the skyrmion will give rise to an addi- 
tional contribution for the Hall voltage beyond that from 
gapped Dirac fermions. 

It seems to be clear that the dissipationless skyrmion 
Hall current will survive on the topological surface state 
in the thermodynamic limit. However, it is difficult 
to guarantee such a phenomenon in the case of non- 
relativistic electrons with the spin-orbit interaction. First 
of all, it should be noted that our effective field theory 
[Eq. ([1])] is valid only in the limit of ^ i > Xso- In 

the Jh limit ps = ''"V'JV'cr^ oc JhS will vanish. 

As a result, we lose terms to describe dynamics of local 
spins in the effective Lagrangian. This results from the 
U(l) projection [Eq. (A6)] of the original SU(2) effec- 
tive theory [Eq. (A5)]. On the other hand, the Jh — 
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limit recovers the Rashba model via an appropriate gauge 
transformation in Eq. (A5). In the strong Hund-coupling 
limit the spin-orbit interaction gives rise to an additional 
internal magnetic field for itinerant electrons, as shown 
in Eq. ([T|) . Inserting the skyrmion configuration into the 
spin-orbit induced gauge field, we can see that the inter- 
nal effective magnetic flux given by the spin-orbit interac- 
tion, / (frlda^izlalpZ/s) - dyizla^g^z/s)], decreases as 

where L is the system size. The effective magnetic 
flux via the gauge field of the skyrmion configuration, 
-pj / (Pr{dxay — dyUx), is also proportional to 1/L^ as 
discussed before. As a result, the spin-orbit interaction 
does not allow the dissipationless skyrmion Hall current 
to survive in the thermodynamic limit when Jh is as- 
sumed to be large in Eq. dH). 

In the Jh — >■ limit the situation is more tricky. As 
well known, the Rashba model {Jh = in Eq. (Al)) 
shows the spin Hall effect [1^, where disorder effects are 
not taken into account. When ferromagnetic interactions 
are turned on, a helical ordered state is expected to ap- 
pear. The question is whether the spin Hall effect sur- 
vives or not in the helical ordered state. When the spin 
Hall effect exists, the dissipationless skyrmion Hall cur- 
rent will be observed in the thermodynamic limit. Unfor- 
tunately, we do not have any definite answer because the 
helical order will change the electron dispersion of the 
Rashba model, which can spoil the spin Hall effect. In 
particular, the ordering wave vector may be incommensu- 
rate generically, making the problem much complicated. 
We leave this interesting problem as a future work. 



III. COMPARISON WITH THE 
LANDAU-LIFSHITZ-GILBERT EQUATION 
APPROACH 

The well known Landau-Lifshitz-Gilbert equation 



dS ^ ^ dS 

is generalized as 



(16) 



ds ^ ^ a-^ ds 
- =.,Bx5--5x- 



S X (js ■ V)S] 



(17) 



in the presence of spin current (lll | . The first term in Eq. 

is the standard precession term with -f ~ gfis/fi > 0, 
where 5' = 2 is the 5-factor and /^_b is the Bohr magne- 
ton. B is an effective magnetic field. The second term 
in Eq. (|16p is the Gilbert damping term, phenomenolog- 
ically introduced, where a is a damping coefficient and 
S is spin. The third term in Eq. pT)) describes the spin- 
transfer torque from the spin current, where js is the 
spin current and a is the Bohr radius. The fourth term 
in Eq. (jl7p represents another torque contribution, per- 
pendicular to the spin-transfer torque, argued to result 



from spin relaxation of conduction electrons. This term 
is called the (3 term due to the coefficient /3. 

One can derive an equation of motion for domain walls, 
vortices, and skyrmions from this generalized Landau- 
Lifshitz-Gilbert equation, resorting to the collective co- 
ordinate method [IJl- Suppose that the right hand side 
of Eq. or Eq. ([T7|) can be derived from an effec- 

tive Hamiltonian -ffe//- Then, one can write down an 
equation of motion for spin dynamics as follows 

dt 5S 
Multiplying S x to both sides, we obtain 



^ dS dS 



dS SJhfi 



where the orthogonality condition S • = is used, 
is called the collective coordinate, representing the core 
position of the skyrmion (domain wall or vortex), where 
i — x,y. Integrating both sides over the two dimensional 
space area A, we reach the following expression for the 
skyrmion dynamics 



4nQ J2 



0=x,v 



'■ dt 



SH, 



eff 



(18) 



where Q is an integer identified with the skyrmion num- 
ber 



^2 5 dS ds 
d rS ■ — — x — — . 



(19) 



eff 



^c\{-idt - T^jco- - Jhti ■ c^Sapcp 



The skyrmion or vortex dynamics has been studied in 
this framework, starting from essentially the same effec- 
tive action as the present model Hamiltonian [Eq. (Al)] 

2m/ 

+S(t>{l - COS0) + —{drnf + V{n)^, 

where n = (sin cos (j), sin sin 0, cos 0) represents a spin 
direction and Vin) introduces spin anisotropy. J is the 
ferromagnetic exchange coupling constant and 4> = dt4>. 

Performing the same unitary transformation as Eq. 
(A3), one obtains the same effective action as Eq. (A5) 
except for the Rashba spin-orbit coupling term 



lat) 



2m 



-JH<T^l^a + 50(1 - cosO) + —{drfif + F(n)|, 

where the Berry gauge field is given hy ar = ■ S = 
—iU^drU. As intensively discussed before, an essential 
effect is given by the spin chirality fluctuation (Berry 
gauge field) to the spin-current minimal coupling term 



HsT - / d'xl ■ ^(1 



cost 
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where the gauge field is represented with two angles. 

Based on this effective action, one can derive an equa- 
tion of motion for the vortex or skyrmion dynamics [27| 



adt^. 



(20) 



G ~ EzS J (P xri- {dxiiy. dyii) — AttSQbz expresses the in- 
ternal magnetic flux given by the skyrmion charge Q. Vs 
represents the velocity associated with the spin current 
of itinerant electrons. U (^) can be regarded as a pin- 
ning potential, and the first term in the right hand side 
describes the associated pinning force, a is the Gilbert 
damping constant. 

Comparing this equation of motion with Eq. (jl7p , the 
extended Landau-Lifshitz-Gilbert equation, we see that 
X S and -^[S x {j, ■ V)S] in Eq. ^ are not 
introduced in this equation of motion. G x dt^ is associ- 
ated with and G x Vg corresponds to j^ijs • V)5'. 

adt^ results from x 

Each term in our Maxwell-Chern-Simons approach 
[Eq. (|4])] has its partner in Eq. ((20)) except for Chern- 
Simons terms proportional to Occ and Qsc- Considering 
the following constituent equations 



y]crj^ = assSe., -(^tjSej, 

^ — ' TT 

cr 



we obtain 



7 T A 2 



„2 H2 4_ _2 r4 „2 2 



from our Maxwell-Chern-Simons approach. On the other 
hand, Eq. pOp gives rise to the following skyrmion mo- 
tion 



{A7rSQ)Hacc/Pei) 
(47rS'Q)2 -t- a2 

{ATTSQ)aiacc/pei) 
(47r5'Q)2 + a2 



Ex 



where the drift velocity of the spin current is given by 
Vg = {(7cc/ PeOExX- We note that the charge current in 
terms of the field contains the contribution of the 
spin current in terms of the electron field c^- In this 
respect the charge conductivity ctcc differs from the actual 
electrical conductivity. Identifying AnSQ with pei(J^g/L^ 
and a with npeiCTss, the above expression becomes 



TTL'^Pel<J^g(Jssa-cc 



-Ex, 



■Ex 



Pel'^ss + TT 7v Pel^^ss 

This result coincides with the first two contributions in 

the MaxwcU-Chern-Simons-equation approach. 

IV. EVALUATION OF CHERN-SIMONS TERMS 

A. Current-current correlation functions 

An important task is to show that the Chern-Simons 
coefficients are non- vanishing. In this section we evaluate 
all current-current correlation functions explicitly. Inte- 
grating over itinerant electrons in £^ of Eq. ([T]), we find 
the local Chern-Simons action [Cos) in Eq- ©• Gener- 
ally, we consider an effective action for gauge fluctuations 



cgauge 
^eff 



dx / dx'i (aM A fx) ) ( nMt(-'-') nr.(x,x') \ ( a,A^') 



where the gauge kernel matrix is given by 

n^t(x,x') = t2^r,{[aJ,,,(x)][a'J,-„(x')]})^-4itA,o(T,{[n^(xK(x)][a'J,v(x')]}) - ^ ^ , p ^ v 

- 16AL(t.{ K(x)p.(x)] [nl (x')p.' (x')] } )^, 
n;^^,(x,x') = t2^r,{[J,^(x)][J,v(x')]})^-4ia,,o(T,{[n^(x)ap,(x)][J,v(x')]})^- (^xf>x,^Oz. 

- 16AL(r. { [nl (x)ap, (x)] [nl (x')aV.' (x')] } 
n^'=,(x,x') = ^2^^,{[aJ,^(x)][J,v(x')]})^-4^^A,o(T,{[7l^(xK(x)][J,v(x')]})^- (xf^x',/xOj.) 

- 16AL(r.{[7i^(x)p,(x)]K(x')aV.'(x')]})^, 



where the summation for ct = ± should be performed. 
""M ~ h'^a'^ap^'^ represents the skyrmion configuration 



(21) 



(22) 



of the spin direction /i, and 

j^^(x) = -0^ (x)[a^i/,^(x)] - [a^V',l(x)]'(/'<T(x), 

p.(x) = Vi(x)V'.(x) (23) 
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denote the current and density with spin a. n*'^(x, x'), 
n^'J,(x, x'), and n^'^,(x, x') are the spin-current- 
spin-current, charge-current-charge-current, and spin- 
current-charge-current correlation functions, respec- 



tively. The subscript c represents "connected." 

Applying the Wick's theorem to the above expression, 
we obtain 



n^t(x,x') = -t" [{9.,a.(x,x')}{a.;^G,.(x',x)} - C(x,x')}g,.(x',x) 

(T=± 

- C(x,x'){5,^a,^G,.(x',x)} + {a,,C(x,x')}{a,,G,.(x',x)} 

-t- AttXsofi'^Xx) ^[G,.(x,x'){a,^G,,(x',x)} - {9,, G„„(x,x')}g,,(x',x) 

<T = ± 

n-(x,x') = n-(x,x'), 

n^t(x,x') = -e Y fT[{a,^G,,(x,x')}{a,. G,,(x',x)} + g,,(x,x')}g,,(x',. 
+ G,,(x,x'){a.,,a,^G,,(x',x)} - {a,, G,,(x,x')}{a,,G,,(x',x)} 

+ AttXson';Xx) Y [G,.(x,x'){a,^G,,(x',x)} - {9,, G,,(x, x')}g,,(x', x) 

(T = ± 

+ 16Xl^n''^{x)nl{x') ^ aGaaix,x')G„„{x\x), 



X •H- X , O I' 



(24) 



cr=± 



where the single particle Green's function is given by 
G,,(x,x') = -(r,{V'.(x)Vt(x')}). (25) 

B. Single particle Green's function 

It is an essential procedure to find the single particle 
Green's function in the presence of the single skyrmion 



r 



due to translational symmetry breaking. An important 
feature that we try to catch is time reversal symmetry 
breaking to allow Hall conductivity, where an internal 
magnetic flux of the skyrmion gives rise to the Hall mo- 
tion of itinerant electrons. This effect will be introduced, 
solving the Schrodinger equation in the single skyrmion 
potential and finding new eigen states. 

We start from the Schrodinger equation with the single 



potential, where the plane wave basis cannot be applied skyrmion potential in the polar coordinate 



{-idt - fJ.r - JHSa)%p„{r, (j), t) 

(( d . 2Xso 2^r \\ If d . 2A,o 2^r 

— ti t(T— — + - ta- 



dr ta + J r\dr ta + J \ + 



1 22 
+ ^(a0 + «(T-j^) }^.,(r,0,t) =0, (26) 



where ^ is the skyrmion core size. In appendix Bl we expand the wave function in a complete basis 
show its derivation. 



Mr,^,t)^ J dS,e-'^'*^a,(n,£;,)*^(r,0;£;,), 



(27) 

where is an energy eigen value and n is a good quan- 
A standard way to solve this kind of equation is to turn number, '^'^{r, (j); E^) is an eigen state described 
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by two conserving quantum numbers, n and E^r, and 
acr{n,E^) is an associated coefficient in the expansion, 
identified with an annihilation operator in the second 
quantization expression. The rotational symmetry of the 
skyrmion potential leads us to the following decomposi- 
tion 



Kir,<P;K)^Ce'''^F^{r;E^), 



(28) 



where n is identified with an angular momentum and 
F^(r;E^) is the corresponding radial wave function. C 
is the normalization constant, determined by 



Unfortunately, we fail to find an analytic expression for 
the radial wave function in the presence of the spin-orbit 
interaction. Instead, wc could obtain the most general 
expression for the radial wave function in the absence of 
the spin-orbit coupling 



drr 



K*{r,cf>;E,)^Ur,c^;E„') 



(29) 



'HeunC(0, V2, - ^f., '-f^ + ^f., 



+ + r^)^r-"HcunC(0, -n, ^2, - ^f., + ^ 



r 



(30) 



where HeunC(o, n, ^/2, - ^£„, + ^£^, -|^) is 

a solution of the Heun's confiuent equation, regarded as 
the generalization of the hypergeometric function when 
an ordinary differential equation contains four singulari- 
ties [l^. £cr = t~^{Ea + iJir + JhSg) IS an effective energy 
level. In appendix B2 we discuss the Heun's confluent 
equation in detail. 

The main feature of this wave function lies in the time 



reversal symmetry breaking. An effective magnetic flux 
due to the skyrmion gives rise to an effective Lorentz 
force for itinerant electrons, resulting in chirality to the 
electron state. One can see this effect from the difference 
between two time reversal states with opposite angular 
momenta. In the asymptotic limit of r — > oo the Heun's 
confluent function becomes 



J 



F^{r)=C[{e + r^'Y-^r^e^^(- 



4.{n+2+ ^/2) 



r 



1 + 



£a 



exp 



4(-n + 2 + y/2) 



as shown in appendix B2. This asymptotic wave function 
reveals that quantum states with negative angular mo- 
menta are suppressed due to the effective Lorentz force, 
imposing the condition that the wave function should not 
diverge in the r — > oo limit. If we consider positive an- 
gular momenta in the above expression, the second term 
allows only n 1, 2 due to the regularity condition, thus 
forcing €'2 to vanish identically. On the other hand, the 
flrst term allows only n = —1,-2 among negative an- 



gular momenta, giving rise to C'l = 0. Considering that 
cither case reaches the same expression, the consistent 
solution is given by only the flrst term with positive an- 
gular momenta. 

It is straightforward to find the single particle Green's 
function in the new basis. Inserting Eq. (|27p with Eq. 
([28| into the Green's function [Eq. ([25]) ] and performing 
the Fourier transformation, we obtain 



J 



G^,(r,r',0-</)',jc^) = |C|^ 



V- Fl^{r-£,)F:i{r'-£,) 



Y,c{iUJ,£c) - £a 



(31) 
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where the self-energy correction 

is introduced, resulting from both elastic impurity scat- 
tering and inelastic interaction effects with gauge fluctu- 
ations, respectively. A way to understand this expression 
is shown in appendix B3. 

C. Longitudinal conductivity 

Inserting the Green's function [Eq. (|3T|) ] into the con- 
ductivity tensor [Eq. ([Ml)] with the following represen- 



A lengthy but straightforward calculation, which is 
given in appendix CI, yields 



where Aa{vij,£a) is the spectral function for itinerant 
electrons. 



respectively. f{r) = -^i^ in S'J{£a;£'a) originates 
from the skyrmion configuration of the spin component. 
We note that only nearest neighbor angular-momentum 
channels are coupled. The cross correlation between the 



tation 



dx = cos (t)dr dd), dy ~ sin (j)dr H dd> 

r r 

(32) 



in the polar coordinate, we are ready to calculate the 
longitudinal "spin conductivity" 



(33) 

I 



(34) 



(35) 



(36) 



spin-current and spin-density, proportional to tXso in Eq. 
([22|) . vanishes for the longitudinal conductivity. 

If we consider the non-interacting limit of 
5ftScr(i'CT), 5So.(i'(t) 0, we can simplify this ex- 



7^ = 
sp 



drr 



dr'r 



2tt 



2tt 



, 5n-(r-,r',0-0',r! + i(5) 

lim . 

n-»-o Q 



16n'Xl\C\'J2 I '^^-j < /^rf^.(-^^)^^^/.(^<x,f.)X(^^.,f;)5^['5"(^-^^)]'' 



1 ^Ti(j{VfjX^) 



and TZ^^{£a-;£'a-) and S^{£cr;£'c^) are composed of radial wave functions, given by 

K{Sy,Sa) - drr([a.F;(r;£.)]F;+i(r;£;) - F;(r;f.)[a.F;+i(r;£;)] 

S:{£,-£'„)= / dTrf{T)F^{r-£,)F:+\r-£'„), 
Jo 
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pression further, resorting to the asymptotic form of procedure is shown in appendix C2, and the final 
the radial wave function in the case of A^o = 0. This analytic expression is given by 



V. 2 ( 71/ I 1 ) 



L^^2F,[2 + n,^V2,3 + n,^{L/0^] 
2{n + 2) 



(37) 



where L is the system size and ^ is the skyrmion core size. The cross longitudinal conductivity is given by 
The hypergeometric function is defined in the Mathemat- 
ica program. 

As discussed before, the electrical conductivity is the 
same as the spin conductivity 

<^ti = <■ (38) 



+ le^^ALlCr 5] cT j d£.j < |_^di.,(-^^)^^^^XK,f,)X(i^.,f;)X^[5;'(£:<.;f;)] 



(39) 



As mentioned before, this correlation function is finite 
because we are considering ferromagnetism. But, this 
cross effect is not crucial for the skyrmion dynamics, just 
modifying the transport coefficient. 



ward calculation given in appendix Dl leads to the fol- 
lowing expression, 



D. Hall conductivity 



The Hall conductivity can be obtained along the same 
strategy as the longitudinal conductivity. A straiglitfor- 

I 



n 

~^nhKo\C\'' Y^a [ dEJ dE', r dv, r dv' J^^^^ ~ ^^''^^ A.{v.,EM.{^'.:Oy^O:,{Ey,E',)S:{E,-E',) 



— oo 
oo 
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where 



poo POO 

/ drF-{r;£,)F-+\r;£'^), V:{£y,£',) = drrF^\r;£,)F-+\r;£'^), 
Jo Jo 

dT([drF:(r- £.m+\r; £',) - F^ir; £,)[5.F;+i(r; £',)]) , 
drr[[drF-{r;£,)]F-+\r;£'^) - F^ ir;£,)[dr>F-+\r; £',)]) , 



drrf{r)F:{r-£,)F:+\r-£',) 



The spin-orbit interaction contributes to the Hall con- 
ductivity. 

This expression reveals that the Hall conductivity does 
not vanish even without the spin-orbit coupling because 
the r integration differs from the r' integration. The 
underlying mechanism is that the chirality is preferred, 
reflected in the Heun's confluent function. As a result, 
the r integral becomes different from the r' integral, orig- 
inating from the angular dependence. 

The finite Hall conductivity becomes clearer if one fo- 
cuses on the asymptotic expression of the radial wave 
function. We simplify Eq. ([Ml further in the non- 
interaction limit, shown in appendix D2. We note that 
the exponential decay of the asymptotic form is expected 



(41) 



r 



to make the contribution from the radial integral finite 
for the Hall conductivity. Since this spin Hall conductiv- 
ity corresponds to af^ = J dx J dyQss{x), the finite Hall 
conductivity means that it is not proportional to as 
discussed before. 

The charge Hall conductivity is the same as the spin 
Hall conductivity 



eH 



(42) 



implying that is a constant, not proportional to 
as discussed in Eq. ([T^ . 

The cross Hall coefficient corresponding to cr^ in Eq. 
dS]) is given by 



= 27Th^\C\' ^ a / < H dv, H di^' J^"^^ _ ^^^"^^ A„{v,.£M.WM 

J J J -OO J -OO l^CT i^o-j 

X Y,[0Z{£.-£X{£a-.£',) - v:{£y,£',)QU£,;£',)} 

11 

+%'kHXso\C\''Y, J d£a j d£'„ J dv, 



dVr, 



d^; ^^!:^'^^^-^,^';^^^ .4.(^.,g.)^.(^;,g;)^0g(g.;£:;)5^(g.;£:;) 



V 



This is also non- vanishing because of the JrS term in the 
effective action, resulting from the different population 
between f and I itinerant electrons, proportional to JhS. 



V. SUMMARY 

We investigated dynamics of skyrmions under spin cur- 
rents driven by electric field in itinerant ferromagnets. 
We developed a novel framework based on the effective 
U(l) gauge theory formulation [Eq. ^ and Eq. ([2])], 
where the Maxwell-Chern-Simons equation [Eq. Q with 
Eq. ©J is the key equation for soliton dynamics. Al- 
though this framework differs from the Landau-Lifshitz- 



Gilbert equation approach, both equations turn out to 
have essentially the same ingredient as it should be. In- 
deed, we recovered the expected result of the Landau- 
Lifshitz-Gilbert equation approach [Eq. (PH])] from the 
Maxwell equation framework [Eq. ([7])]. 

An important improvement beyond the previous study 
lies in the mutual feedback effect for both skyrmion dy- 
namics and electron motion, where the internal magnetic 
fiux of the skyrmion gives rise to the Hall motion to 
itinerant electrons or vice versa. This physics can be 
described by the Chcrn-Simons terms [Eq. ^] in the 
skyrmion moving frame. As a result, we revealed that 
the skyrmion motion follows not only the electric field 
but also its transverse direction [Eq. PH)) ]. An interest- 
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ing observation is that even if an insulating state is con- 
sidered, the electric field will induce the dissipationless 
skyrmion Hall current due to the Chern- Simons terms 
[Eq. ([T^ ]. In particular, we predict that the dissipa- 
tionless skyrmion Hall current will survive even in the 
thermodynamic limit as far as the spin-orbit interaction 
is introduced, expected to realize in the surface state of 
three dimensional topological insulators when magnetic 
impurities are deposited. 

We evaluated the Hall conductivity in the presence of 
the single skyrmion potential, showing that it is non- 
vanishing indeed although it will vanish in the thermody- 
namic limit because we are considering only one skyrmion 
without the spin-orbit interaction. A careful treatment 
is required because we should construct new eigen basis 
in the single skyrmion system [Section IV-B]. Based on 
this construction, we calculated all kinds of correlation 
functions explicitly, and found the general expression for 
the HaU conductivity [Eq. (gD]) with Eq. (|IT|) ]. 
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search under Grant No. 19052006 from the MEXT of 
Japan and No. 21740275 from Japan Society of the Pro- 
motion of Science. 



Appendix A: Mapping from the double exchange 
model to an effective U(l) gauge theory in the 
strong coupling limit 

We start from the ferromagnetic Kondo lattice model 
or the double exchange model with the Rashba spin-orbit 
coupling 




i ij 

■i 

-iKo^<fi{c\+^o.c^p - cL^+a/g), (Al) 



where Cia represents the conduction electron field and Si 
expresses the localized spin, /i is an electron chemical 
potential and t is the wave-function overlap integral for 
conduction electrons. Jh is the exchange coupling con- 
stant, set to be positive. A^o is the Rashba spin-orbit 
coupling constant, which breaks the inversion symme- 
try, realized on the surface or in the inversion-symmetry 
breaking material. Sb is the single-spin Berry phase term 
in the spin coherent-state representation, given by 



r0 

Se^iS dT^dr(j),{l-cose,), (A2) 



where the spin field is expressed by two angles. Si — 
(sin Oi cos 4>i , sin 6i sin (pi , cos 9i ) . 

If we consider the system such that dynamics of local- 
ized spins is much slower than that of itinerant electrons, 
spins of conduction electrons follow localized spins. This 
situation allows us to take the strong coupling approach 



— JHSi-aaf) — —JHSUiafO-^s^iSfi' 

Ipia = Uj^pC^f}, (A3) 



where the unitary matrix field Ui = \ * *i | con- 

sists of a bosonic spinon field Zi^ and the alignment 
of the itinerant spin to the localized spin introduces 
a electron field ipia. The spinon field can be repre- 
sented in the following way, z^-j- = e^*"?" cos ^ s-nd 
Zii = sin^, which gives the self-consistent expres- 
sion Si = ^zj^^apZip. 

Representing Eq. (Al) for Cig- and Si in terms of V'm 
and Uia/3 [Eq. (A3)], and taking the continuum approx- 
imation, we reach the following expression 
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+ti;Ud,U^ + [U^drU]uA (drU + U[{drU^)U]) - t^r (^I ) 

\ / a-y V / 7/3 

(A4) 

I 

Introducing the Berry gauge field as the spin connee- as foUows 
tion, A'^p = —i[{dvU^U]aj3, we can rewrite the above 

I 

+t(d.4 + ii^i[Alp + {Koiml^a^^sUsp]) (9r^0 - ^[A^p, + {Kolt)Ul^alsUs^]i,^) 

I 



Remember that this effective field theory is just the 
change of variables in the microscopic model Eq. (Al). 

Because this SU(2) gauge theory formulation is quite 
complicated, we perform the U(l) approximation. One 



can understand this procedure as the staggered-flux 
ansatz in the SU(2) slave-boson theory, where the SU(2) 
gauge symmetry is reduced to the U(l) symmetry jl9j . 
Then, we find an effective U(l) gauge theory [2^ 



Z = j Di',Dz^Da^6{\z,\^ - l),5(c),a,)e-^^"^'' '^-J'^''-^, £ ^ C^, + 

= Psztidr - iar)za + tps\[dr - la^ - i{Xso I t)zl^al^iiZp\zc\^ , (A6) 



where an electromagnetic vector potential Ar is intro- 
duced. It is interesting to see that the internal gauge field 
Or couples to the spin current of the electron field while 
the electromagnetic field Aj- does to the charge current, 
/i^ = ^ + 2t + 2A^o/i is an effective chemical potential 
for itinerant electrons. The spinon part is reduced to the 
CP^ representation of the ferromagnetic 0(3) nonlinear 
cr model, where the time derivative term is added ex- 
plicitly. This time derivative term is expected to appear 
from quantum corrections, i.e., the self-energy correction 
to the spinon dynamics. In the antiferromagnetic case 
the second order time derivative term can be applied to. 
The Rashba spin-orbit coupling gives rise to an interac- 



tion term between the spin current and spin, quenching 
the spin direction to the momentum or current direction. 

Considering that the Berry phase term is associated 
with a background potential, we arc allowed to take into 
account the saddle-point configuration of the gauge field 
for the Berry phase term 

= -^[4(5m^<t) - (5m4)2t] = cos 61. 

Then, the Berry phase term can be written as follows 

SB=2iS dr J (frar, (A7) 
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where iS dr J (Prdr4> = is used for the skyrmion configuration. 



J 



Appendix B: Single particle Green's function in the single skyrmion potential 
1. Schrodinger equation with the single skyrmion potential 

We start from the Schrodinger equation with the single skyrmion potential 

{—idt — /ir — JhSct — i(7ar)lpcy — t{dr — ICJa-r — lA^Y'^-'a = 0, (Bl) 

where the spin-orbit interaction is neglected. Introducing the polar coordinate associated with the symmetry of the 
skyrmion potential, 

-t(fdr + 0^ - i(Ja{r)^^ ■ {rdr + 4>y ^ icTa{r)4>^ = -t{dl + ^ + ^ - 2iaa{r)^ - a^r)'^ , (B2) 
where the gauge potential a{r)4> = ax{x, y)x + ay{x, y)y in the polar coordinate is given by 



a{r) = ^al + al^^^ (B3) 
for the single skyrmion solution, we rewrite the above Schrodinger equation as follows 

{-idt- iJLr- JHS(T)i;„{r,(t),t)-t\^dl + ^ + --^{d^ + i(T-^^^^ }^/'^(r, </>, i) = 0. (B4) 

2. Heun's confluent equation 

Inserting Eq. ([27]) with Eq. ([28]) into Eq. (B4), we obtain the eigen value problem for the radial part 

dlF^ir) + ^F:{r) -^{n + a-^^)" F^^{r) + r^E, + + JHSa)F^{r) = 0. (B5) 

Performing the change of variables = and introducing F^{t) = t'P{t — lYY^{t) with constants p and g, this 
equation can be written as follows 

where Ea = t^^{Ea + jUr + JhSct) is an effective energy level. Surprisingly, this equation is known to be the Heun's 
confluent equation, and its solution is well understood as follows 

i;"(t) = HeunC(0, n,V2,^^ - ^£,, + ^f., t) . (B7) 

A general expression of the Heun's confluent equation [28j is given by 

(fY{z) [-az^ + {a- P --f -2)z + P + l] dY{z) [{-^ - 7 - 2)a - 25]z + {P + l)z + (-7 - l)/3 - 7 - 2r; 
dz2 z{z- 1) dz 2z{z - 1) 



Y{z) 



satisfying two boundary conditions such as 

(B9) 



Y{z = 0) = 1, 

dY{z) (-a + 7 + l)/3 + 7-a + 2?7 



dz 



z=o 2(/3+l) 
A general solution of this equation is known to be the Heun's confluent function 

y(z) = HcunC(a,/3,7,(5,?7;z), (BIO) 
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where p and q are determined by 

1 

= 9(<Z-1) = 4, (Bll) 

and all other coefBcicnts arc given by p and q in the following way 

a = 0, /3 = 2p, 7 = 2g-l, 5 = p+--^£:., = + (B12) 

Two independent quantum numbers appear to be an angular momentum n and an energy eigen value £cr, where n is 
an integer while £„ turns out to be continuous. 

In order to understand time reversal symmetry breaking in the Heun's confluent function, it is valuable to find its 
asymptotic form in the r — >■ oo limit. The corresponding Schrodinger equation is given by 

2p + 2g+l g^y„(^^ _ l^f.C(t) « 0, (B13) 

where dominant terms are selected by the t ^ —oo limit. It is straightforward to solve this differential equation. 
Introducing p = n/2 and q = (1 + y/2)/2 into the solution, we obtain 

y"(r c») cx expf — ^rA . (B14) 

' V 4(71 + 2 + ^2) / 

3. A Green's function in a skyrmion background 

Inserting Eq. ([?7)) with Eq. ([25]) into the Green's function, we obtain 

C(r,r',0-0',t-O = -|C|2 / d£,e-'^^-(*-*')^e'"(^-^V;(r;f,)F;(r';f.)(a<,(n,£:<,)at(n,£<,)). (B15) 

n 

We note that the radial coordinate cannot be r — r' due to translational symmetry breaking. 
Performing the Fourier transformation for time, we obtain 

G<..(r,r',0-0',.. + z<5) = |Cp / ^ e^^^-^V; (r; £.)F; (r'; ( + ^^/%^) (B16) 

in the real frequency, where f{£a) is the Fermi-Dirac distribution function, while given by 

GUr,r',^-^',^UJ) = \C\' f dg^y^^.ni^-^')ElM^}I^^I^ (B17) 

in the Matsubara frequency. 

We introduce the self-energy correction, resulting from both elastic impurity scattering and inelastic interaction 
effects with gauge fluctuations 

= Er^'li^^f^) + (B18) 
Then, the most general expression for the single particle Green's function is given by 

J ^ lUJ - L,a{iUJ,C.a) - Ca 

Appendix C: Polarization functions for longitudinal conductivity 
1. Formal expressions 

For convenience, we decompose into four contributions; 

n-(r,/,0-(/.',*r!)= n-W(r,r',0-(/.',zf}). (CI) 

K=A,B,C,D 
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-F, 



The first part is the conventional particle-hole channel in the presence of the single skyrmion potential, 
n^;(^)(^'''>-0',«^^) = -i'|C|'*cos0cos0'^^e'("-"''(^-^') / d£„ f d£'^ 

n n' J J 

{[drF-{r;£^)]F:{r';£,)[dr,Ff{r';£'J]F-\r;£'^)~[drF-{r;£,)^^^^^ 
-F^{r;£,)F:{r';£,)[dr'F:\r';£',)][drF'^\r;£',)] + F:,^ir;£,^^^^^ 

r dv, r di^;iMi:i^lx(i^.,fa)X(z^;,f;), (C2) 

where the spectral function for itinerant electrons is given by Eq. psp . In the ideal non- interacting limit 
3fiEo-(z/cr), 3Scr(t'(T) — we obtain 

r dv. r dv'^i^^i^^AAvMA.w.K) ^ .y^L'^^^il 

J-oo J-oo lil-{Va-v'^) in-{£„-£'^) 

nothing but the particle-hole polarization function in the Fermi gas. 
Other pieces are given by 

^ri^)(r,r^</)-0^^^7) = ^2|c|4!i]^cos(/.'^^^(n-n')e^("-"')(*-^') j d£„ f < 

^ n n' J J 

{F^ir; £.)F^ir'; £.)[dr,F:' (/; £',)]Ff (r; £',) - F^{r; £,)[dr>F^{r';£,)]F^' (r'; £'„)F^^' {r; £',) 
:{r; £^)F-{r'; £,)F-' (/; £'„)[drF-' {r- £'^)] + F;(r; £^)Fl\r'- £^)Ff {r' ■,£'J[drF-' (r; £'^)] } 

d'^a r dv'„ .(j'^^l" ^^^V (C4) 

n-(^)(r,r',0-</)',zl7) = -<2|C|4cos0^^^*(n-n')e'("-"')(^-^') / d£, f < 

{ [drF-{r; £.)]F-{r'- £^)F-' (r'; £'„)F:' (r; f^) - [a.F;(r; (/; (r; £'„) 

-F:{r- £.)F^ir'; (r'; £',)]Ff (r; 4) + F;(r; £,)[dr>F^ir'; £,)]Ff (r'; (r; £',)} 

r dv, r dv',I^^^^-^^A.{v,.£,)A.{v',,£'„), (C5) 

and 

n n' J J 

[F:{r- £.)F:{r'- £,)F:' (/; f^jF^' (r; f^) - F;(r; (/; £X' (r; f^) 

-F;(r; £.)F;(r'; (r'; (r; f^) + F;(r; (r'; (r; £',)} 

r dv^ r d,.'^^^^]^i^\AA'^a,£.)AAK,c)- (C6) 

J-oo J-oo - {^<y - ^a) 

Substituting the above expressions into Eq. ([33]), it is found that the contributions from IlJt^^'' (r, r', — (f)',in) 
and Uf^'-^\r,r',(l)- (t)',in) +UfJ-'^\r,r',(j)- (j)',in) all vanish. Therefore, we obtain Eq. ((3i)) . 

2. Asymptotic form 

The conductivity is simplified as follows in the non-interacting limit without the spin-orbit interaction 
a%^2nh'\Ctj d£^[-^J^)^^^J£[ni{£^-£^)f. (C7) 



18 



Resorting to the asymptotic expression, we see 

drF-{r ^ oo;£,) = ((1 + V2)-^ + - - ^^r)F;(r ^ oo;£„). (C8) 

V + H r 2{n + 2 + V2) ' 

Then, we obtain 

nii^EM = - / drr(- + ^ r) f;-(r; g.)^^ ^(r; g.) 

Jo 2(n + 2 + \/2)(n + 3 + V2) / 

^o- \,.-?. , 2\1+a/2 2n + 1 / 2n + 5 + 2-\/2 



2(n + 2 + V2)(n + 3 + V2) / ' V 4(n + 2 + V2)(n + 3 + V2) 

(C9) 

Inserting this expression into Eq. (C7), we reach the following expression in the T — > limit 



2(n + 2 + \/2)(n + 3 + V2) 
' H "4(„ + 2 + y2)(n + 3 + v^) 



n "'0 

L^2^fi[2 + n,-V2,3 + n,- (L/e)^] j^ ^^^^^^ 



2(n + 2) 

where L/^ is the ratio between the system size and the skyrmion core size and 2^1 is the hypergeometric function. 

Appendix D: Polarization functions for Hall conductivity 
1. Formal expressions 

It is convenient to decompose 11*^ into four contributions; 

Klir,r',^^cj,',zn)= E n-(^)(r,r',</>-0',zl7). (Dl) 

K=A,B,C,D 

The four polarization functions are given by 

nj;(^)(r,r',0-</)',iO) = -<2|C|4cos0sin0'EE^'^"~"'^^*~*'^ / / 

n n' 

{ [drF^\r; £.)]F^ir'; S^M^F^^' (r'; E'^W^' (r; - [5,F;(r; £,)] [dr>F^ir'; £,)]F^' (r'; £',)F^' (r; £',) 
'(r; £.)F:ir'; £.)[dr,Ff (r'; f^)] [^.F^' (r; 4)] + ^^(r; f [a.,F;(r'; f,)]^/ (r'; f^) [a.F^ (r; f ;)] } 

n-(^)(r,r',0-(/.',*17) = i2|C|4^sin0'EE*('^-"V^"""'^^^"^'^ / / < 
{f:{t- £.)F:{t'- £.)[dr.F:' (/; £'^)]F:' (r; f^) - F;(r; f.)[9.,F;(r'; (r'; (r; f^) 

-F; (r ; (r' ; (r' ; ) (r ; f ^ )] + FJ' (r ; )F; (r' ; )F;' (r' ; £'„ ) [a.FJ'' (r ; )] } 

/°° dv, r dv',I^^^^^^^^Aa{v,,£a)A.{v',,£',), (D3) 



-F 
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^ n n' J J 

{[drF^ir; E,m{r'-E,)F^' (r'; E^' (r; E',) - [drF^ir- E„m{r' ■E^)F:' (/; E^' (r; E',) 
-F:{r- E.)F:{r';E,)[dr>F:'{r'-E',m' (r; 4) + F;(r; f ,)[a,.F;(r'; (/; 4)F;'(r; E',)] 

dv^ r dv'^^^^^^fl^AAv^,E,)A„{v'„,E'^l (D4) 



and 



n n' 

{F;(r; £.)F;(r'; (/; (r; f^) - F;(r; (/; (r; f^) 

-F;(r; f.)F;(r'; f.)^^;' {r';E'^)F-' (r; f^) + F;(r; f.)^^;' (r'; (r; f^)} 



oo 



diy. I du',^p^-J^^A.{v.,E,)A„K,E',). (D5) 



It is found that li-^xy ^^ {r, r' , 0—0', iVL) and n^t/^^-* (r, r', 4>~4>' , iVL) do not contribute to uf^, while Hx^y ^\r, r' , 0—0', if!) 
and nrj^)(r,r',0-0',iri) ffive finite contributions. 

2. Asymptotic forms 

In the non-interacting limit the Hall conductivity is given by 

atH^27rh^\C\^l dE^j dE'J-^^^—^Y.{^''^^^'^-^^X^^^^^ (D6) 

where the spin-orbit interaction is not introduced. 

Inserting the asymptotic expression of the radial wave function as performed in appendix C2, we obtain 

V 1.4(n + 2+v/2) 4(n + 3 + %/2)J / 

°° -1 & a 



Jo 



2{n + 2 + V2) 2(n + 3 + \/2) 

X (C2 + ^2)i+y2^,2„+i /_r £^ £^ ^jr^y (D7) 

^ ' V U(n + 2 + V2) 4(n + 3 + V2)J / 

It is difficult to perform further simplification analytically for the general case of n. However, it is clear that this 
expression does not vanish due to the factor of r in the r in tegration. 
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